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We have measured the X-ray spectral properties of a complete flux-limited 
sample of bright AGNs from HEAO-1 all-sky catalogs to investigate their statis- 
tics and provide greater constraints on the bright-end of the hard X-ray lumi- 
nosity function (HXLF) of AGNs and the AGN population synthesis model of 
the X-ray background. Spectral studies using data from ASCA, XMM-Newton 
and/or Beppo-SAX observations have been made for almost all AGNs in this 
sample. 

The spectral measurements enable us to construct the neutral absorb- 
ing column density (log A^h) distribution and separate HXLFs for absorbed 
(log Nn[cm~'^] > 21.5) and unabsorbed AGNs in the local universe. Our re- 
sults show evidence for a difference in the shapes of HXLFs of absorbed and 
unabsorbed AGNs in that absorbed AGN HXLF drops more rapidly at higher 
luminosities than that of unabsorbed AGNs, which is similar to that previously 
reported. In the Lx — A^h plot, we found no AGN in the high-luminosity high- 
intrinsic absorption regime (log Lx[erg/s] > 44.5, log A^H[cm^^] > 2L5) in our 
sample, where we expect ~ 5 AGNs if we assume that absorbed and unabsorbed 
having identical AGN HXLF shapes. We also find that the observed fiux with 
ASCA or XMM-Newton is smaller than that with HEAO-1 by a factor of 0.29 
on average, which is expected for re-observation of sources with a factor ~ 2.5 
variability amplitude scale. 

Subject headings: galaxies: active — quasars: general — surveys — X-rays 



1. INTRODUCTION 

Investigating the X-ray luminosity function (XLF) of AGNs in different redshifts is very 
important in understanding the growth of super massive black holes (SMBH) at the centers 
of galaxies over the cosmological time scale, which is directly linked to the accretion history 
of the universe. Recent progress on X-ray surveys has enabled us to trace the evolution 
of the XLF over a wide redshift and luminosity ranges. In particular, in the soft X-ray 
(0.5 — 2.0 keV) band, a combination of extensive ROSAT surveys (Miyaji et al. 2000b) has 
been extended by recent Chandra and XMM-Newton surveys (Hasinger et al. 2005). With 
this combination, consisting of ~ 1000 type 1 AGNs, they traced the evolution of the soft 
X-ray luminosity function (SXLF) with a high precision. In particular, they clearly traced 
the so-called anti-hierarchical AGN evolution (some call 'down-sizing'), where the number 
density peak of lower luminosity AGNs comes much later in the history of the universe 
{z ~ 1) than that of higher luminosity ones (^ > 2). A fundamental limitation of the soft 



-3- 



X-ray sample, however, is that it is not sensitive to type 2 AGNs, in which AGN activities 
are heavily obscured by photoelectric absorption by intervening gas, except those at very 
high redshifts. 

Hard X-ray (2 <£■<? — 10 keV) surveys are sensitive also to type 2 AGNs, up to a 
column density of A^h ~ 10^^ cm~^ (Compton-thin), giving a much more direct measure of the 
accretion history onto SMBH. Several groups constructed hard X-ray luminosity functions 
(HXLFs) of AGNs in this band (Ueda et al. 2003; Barger et al. 2005; Silverman et al. 2005; 
La Franca et al. 2005). Hard X-ray surveys are limited in the number of AGNs as well as 
completeness, partially due to lower sensitivity of the instruments to harder X-rays and also 
due to the difficulty of identifying optically faint type 2 AGNs. In spite of these difficulties, 
these HXLF studies came to a point showing the anti-hierarchical evolution trend. 

The XLF and X-ray spectra of the constituent AGNs are key elements of the "AGN 
population synthesis modeling" , initially aimed at explaining the Cosmic X-ray Background 
(CXRB) (e.g. Madau et al. 1994; Comastri et al. 1995; GiUi et al. 2001), where the contri- 
butions of cosmologically evolving populations of AGNs with different absorption columns 
(Ah) are considered to synthesize the CXRB spectrum, source number counts and the dis- 
tribution of column densities in X-ray surveys with various depths, areas and energy bands. 
In order to test/constrain these models, a complete samples with well defined flux limits 
as a function of survey area is essential. Investigating the spectral properties such com- 
plete samples models will be particularly strong observational constraints by providing, e.g. 
the distribution of A^h- Ueda et al. (2003) (hcaraftcr U03) constructed an HXLF in the 
2—10 keV luminosity range of logLx (erg s^^) 41.5 — 46.5 as a function of redshift up to 
3, using an extensive set of highly complete samples (> 96% total) of 247 AGNs from Hard 
X-ray surveys from HEAO-1, ASCA, CHANDRA. U03 also integrated the X-ray spectral 
(or hardness) information, with the intrinsic absorption column density as a main spectral 
parameter into the HXLF to construct a population synthesis model. One of the important 
findings of U03's analysis is that the ratio of obscured to unobscured AGNs decreases with 
luminosity, as (less quantitatively) suggested by Lawrence & Elvis (1982) and Miyaji et al. 
(2000a). This has been verified in the nearby AGNs by Sazonov (2004) from a sample of 
95 AGNs detected in the 3 — 20 keV band during the RXTE slew survey (Revnivtsev et al. 
2004). 

In order to investigate the nature and evolution of the X-ray emission of AGNs, sta- 
tistical characteristics of the X-ray sources in the bright-end sample (mostly consisting of 
present-day population) from large area surveys is important to compare with the fainter 
(high-^) population from deep surveys. The sample has to be complete with well-defined 
criteria, so that we can adequately account for selection biases. The spectral properties of all 
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sources of the sample should be analyzed with good spectral resolution data (for example, 
about 1000 cts with 7 ks of XMM-Newton / 1000 cts with 20 ks of ASCA for 5e-12 cgs) 
to constrain F to d=0.03 and logiVn to ±0.3. We also note that there are many absorbed 
AGNs, the accurate measurements of A^h enables us to calculate de-absorbed luminosity, 
which is a more direct indicator of the intrinsic AGN power and thus the mass accretion 
rate. One such study has been made by Schartel et al. (1997) for Piccinotti et al. (1982) 
sample from HEAO-1 A2 using ROSAT ojid EXOSAT as the source of spectroscopic infor- 
mation. A recent study of RXTE Slew Survey catalog (Revnivtsev et al. 2004) utihzed by 
Sazonov (2004) has absorption estimates based on hardness ratio in two RXTE bands, but 
this method is only sensitive to neutral absorbing column densities of log Ah (cm~^) > 22. 
In deeper regimes, Mainieri et al. (2005); Mateos et al. (2005a, b) made extensive spectral 
analysis on the X-ray sources detected in a XMM-Newton serendipitous and deep Lockman 
Hole surveys, sampling intermediate to high redshifts. In view of these, we have compiled 
the results of high-quality intermediate-resolution spectroscopic analysis for a complete hard 
X-ray fiux-limited sample of AGNs defined from HEAO-1 all-sky survey catalogs. Detailed 
spectral information of almost all of the AGNs in this sample have been obtained by making 
spectral analysis and/or by referring to the spectral fit results from literature using data 
from XMM-Newton, ASCA, and Beppo-SAX. The results of our prehminary analysis has 
been integrated into the global HXLF analysis over cosmological timescales by U03. In this 
paper, we present full results of our spectral analysis with necessary refinements and updates. 

The scope of this paper is as follows. In Sect. 2, we explain the construction of our 
sample. In Sect. 3, we explain the spectral analysis of the XMM-Newton and ASCA data. 
The statistical properties of the spectral results, including the distributions of spectral pa- 
rameters are presented in Sect. 4. In Sect. 5, we construct the local HXLFs separately for 
absorbed and unabsorbed AGNs from our sample and show the decrease of the fraction of 
absorbed AGNs towards high luminosities. The overall results arc discussed in Sect. 6. In 
Appendix A, we explain our first-order correction to count rate and the effective survey area 
function for a bias due to confusion noise to the HEAO-1 Al sample. 

Throughout this paper, we use cosmological parameters (Hq, ^a)= (70 hygkm s~^ 
Mpc-^ 0.3, 0.7). When the units are omitted, is measured in erg s and is in cm . 
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2. SAMPLE SELECTION 
2.1. Sample 1 

As a part of our sample (Sample 1), we have used emission-line AGNs from Piccinotti 
et al. (1982), with updated identifications in the on-line catalog (A2PIC) provided by the 
High Energy Astrophysics Science Archive Research Center (HEASARC)^. We only use 
AGNs brighter than 1.25 (R15c/s) in the first scan, which defines their complete fiux-limited 
source list. Only one of the 66 sources above this count rate remains unidentified. The limit 
corresponds to ~ 2.7 x 10^^^ erg cm~^ s^^ in the 2 — 10 keV band for a power-law spectrum 
with photon index 1.65. In the Piccinotti et al. (1982)'s catalog, the regions between -20° 
and -1-20° in Galactic latitude and the Large Magellanic Cloud (LMC) have been excluded to 
minimize contaminations from Galactic sources. Also two AGNs which are within the region 
of the sky covered by Sample 2 (see below) have been excluded from Sample 1 (NGC 4151, 
NGC 5548). As a result of this exclusion, our Sample 1 consists of 28 AGNs with a survey area 
of 2.2 X lO'' deg2. Malizia et al. (2002) pointed out that the z = 0.167 AGN which was listed 
in the on-line A2P1C catalog was a misidentification of the X-ray source H 0917-074. Their 
Beppo-SAX/NFl observation found that the Seyfert 2 galaxy MCG -1-24-12 at ^ = 0.0198, 
which lies within the original HEAO-1/ A2 error box of H 0917-074, dominated the X-ray 
emission in the FOV and had a much larger flux than the z = 0.167 AGN. Therefore they 
concluded that MCG -1-24-12 was the true counterpart of the HEAO-1 A2 source. Thus we 
take MCG -1-24-12 as the identification of H 0917-074. 



2.2. Sample 2 

In order to increase the number of objects in the sample, we have also utilized a deeper 
catalog of X-ray sources detected with the HEAO-1 Al experiment by Wood et al. (1984). 
Identifications of these sources with the help of HEAO-1 A3 experiment (modulation colli- 
mator) have been integrated as the "MC-LASS catalog" by R. Remillard. Both catalogs are 
available on-hne from HEASARC (named Al and A3 respectively). Based on these catalogs, 
Grossan (1992) defined a fiux-limited subsample of 96 AGNs (the Large Area Sky Survey/ 
Modulation Collimator identified sample of AGN or the LMA sample), in which he selected 
those with >0.0036 LASS cts s^^ cm^^, where the unit is defined in Wood et al. (1984). 
Grossan (1992) defined this limit to make a complete sample in |6| > 20°. 



^http: / /heasarc. gsfc.nasa.gov/ 
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Since it is not practical to make intermediate-resolution X-ray spectroscopic foUowup 
for all of the 96 AGNs, we have decided to limit the region for constructing a complete 
fiux-limited sample. For this purpose, we selected a 55 degree radius region from the North 
Ecliptic Pole. This region was chosen to include the maximum number of AGNs which 
had already been observed or planned to be observed by ASCA or XMM-Newton, i.e., to 
minimize the number of new observations, at the time of our initial ASCA AO-7 proposal. 
The defined area covers 25 % of the | 6 |> 20° sky (5.1x10^ deg^). Among the LMA AGNs 
within this region, we have excluded 3C 351 from our sample because it has been found to 
be confused with a BL-LAC object with a comparable flux from our ASCA observation. As 
a result. Sample 2 contains 21 AGNs brighter than the HEAO-1 Al count rate greater than 
0.0036 LASS c/s, which correspond to ~ 1.9 x 10"^^ erg cm'^ s^^ in the 2-10 keV band. 
The overall identification completeness at this limit in our defined area is 85%. However, 
Grossan (1992) argues, based on the comparison with Einstein and i^OS'^T data, that most 
of the unidentified sources are Galactic stars (active coronac) or BL-Lac objects. 

Due to large beams and low count rate threshold, the HEAO-1 Al count rate is subject 
to a bias due to confusion noise. We have made a first-order correction to this effect, as 
described in Appendix A. 

2.3. X-ray Spectroscopic Data and Sample Summary 

The observational data summary of each sources are listed in Tables 1 & 2. All of the 
AGNs in samples 1 and 2 are classified optically as emission-line Seyfert galaxies or quasi 
stellar objects (QSOs). Those classified as BL Lac objects are excluded from our sample. 

In order to obtain good quahty spectra of all the AGNs in our samples, we observed the 
AGN which have not been previously observed nor scheduled with ASCA, XMM-Newton, nor 
Beppo-SAX. We have made our own observations of five AGNs (Kaz 102, 3C 351, MKN 885, 
H 1318+692, H 1320+551) with ASCA (A07) and 2 AGNs (MKN 464 and H 0917+074) 
with XMM-Newton (AOl) to almost complete the spectroscopic study of the sample. Our 
ASCA observation of 3C 351 and XMM-Newton data on H 0917+074 have not been used 
because of the reasons shown above. We have used the spectral analysis result of MCG -1- 
24-12 with the Beppo-SAX ohseTvation by Malizia et al. (2002) instead of our XMM-Newton 
observation of H 0917+074. 

Figures 1 and 2 show the spatial distribution and the relation between X-ray observed 
luminosity and redshift of all the samples. Putting both samples together, we have analyzed 
21 sources of XMM-Newton observation data (13 sources of Sample 1, 8 sources of Sample 2) 
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and 14 sources of ASCA observation data (3 sources of Sample 1, 11 sources of Sample 2). 
For the other sources, we have used spectral analysis results in literature. Among them, 
5 sources were observed with XMM-Newton, 6 with ASCA, and 2 with Beppo-SAX. Thus 
we get the spectral information on all but one (H 1537+339) AGN in the sample of 49 
AGNs defined above. There is no ASCA, XMM-Newton or Beppo-SAX observation data for 
H 1537+339. Remillard et al. (1993) showed an optical spectrum of the AGN identified with 
this X-ray source. It had broad-hnes and a UV excess, which were typical of a type 1 AGN, 
and showed no sign of broad absorption hnes (BALs). Thus we have assumed that it has an 
unabsorbed power-law spectrum (F = 1.8), which is typical of a type-1 AGN/QSO. 



3. DATA REDUCTION AND SPECTRAL ANALYSIS 

3.1. Overall Procedure 

We have reduced the XMM-Newton data using the Scientific Analysis System (SAS) 
ver 5.4.1. or later. We have compared the results of the reduction/analysis using SAS ver 
6.0 and "Current Calibration Files" (CCF) updated as of December 2004 with those using 
SAS 5.4.1. for a few selected sources. The change of the results was insignificant for our 
purposes. About a half of the sources with available XMM-Newton data have such high 
count rates that CCD pile up affects the spectral analysis. In these cases, we extracted the 
source spectrum from an annular region with inner and outer radii of about 12" and 60" 
respectively (exact radii vary) to avoid the central peak. A background spectrum was taken 
from an annular region around the source with the inner radius of 1', when the observation 
was made with the full window or large window mode. Small window mode observations 
are usually made for bright sources and the background are negligible for these observations. 
Since the XMM-Newton light curve of NGC 5548 had a fiare in the last 1/4 of the observation, 
we did not include this time interval into spectral analysis. In the spectral fitting, we used 
pulse-height channels corresponding to 0.4 — 10 keV and 0.3 — 10 keV for MOS and FN 
detectors respectively. 

We have used FTOOLS 5.0 and asca-ANL ver 2.00 to make data reduction and ex- 
traction of events for ASCA observation data. The SIS and CIS pulse-height spectra have 
been created from events screened with the standard screening criteria. Source spectra has 
been extracted from a 3' radius circle and background spectra has been accumulated from an 
off-source area of the same observation for SIS spectra. For the GIS spectra, we have used 
a 7' radius region as a source and the backgrounds from 8' and 11' annulus region centered 
on each source of the same observation. The ASCA light curve of MKN 1152 had a flare 
in the middle and we have excluded the data during the flare. In the spectral flts, we have 
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used pulse-height channels corresponding to 0.7 — 10 keV and 1.0 — 10 keV for SIS and GIS 
detectors respectively. 

For both XMM-Newton and ASCA data, we have used XSPEC ver 11.0 or later for 
the spectral analysis. As our first step, we have fitted the reduced pulse-height spectra 
with a single power- law and intrinsic absorption model (XSPEC model "zwabs"). Galactic 
absorption ("wabs") with a fixed column density (Dickey & Lockman 1990) has been always 
included. This first step analysis have been made in the 2.0 — 10 kcV range to avoid the 
effects of soft excess. We have also excluded the 5.5 — 8.0 keV range when there was a 
significant residual caused by absorption edge or 6.4 keV Fe Ka line in the spectrum. The 
second step has been to add additional components, including a soft excess, an absorption 
edge, a reflection component, or a 6.4 keV Fe Ka line to the spectral model and fitted to 
the full energy range whenever needed to verify that the deviation from the simple best-fit 
absorbed power-law does not have major effects on the spectral statistic and luminosity 
function analysis which arc made later in this paper. We have also tried to include partial 
covering and/or ionized absorbers for some sources as described in Sect. 3.2 . 

The basic results of the simple absorbed power-law fit parameters (photon index F and 
intrinsic neutral absorber's column density A^h) are summarized in Tables 3 & 4. These 
numbers arc from our first-step analysis or similar simple absorbed power-law fits from 
literature except for those detailed in Sect. 3.2. 

Piccinotti et al. (1982), used the conversion factor of one R15 count s~^ ~ 2.17 x 10~^^ 
erg cm~^ s""*^ in the 2 — 10 keV energy band for a power-law spectrum with photon index 
1.65, to derive the fiux for most of their AGNs. Since wc have spectral information for each 
object, wc have re-calculated the fiuxes of Sample 1 objects from the R15 count rate using the 
actual HEAO-1 A2 response function and our spectral results. The corrected observed flux is 
listed as the HEAO-1 flux in Tables 3& 4. An underlying assumption is that F, intrinsic A^h 
and other spectral parameters (except the global normalization) of an AGN do not change 
with the variation of its luminosity. For Sample 2, Grossan (1992) used conversion factor 
5.29 X 10^^ erg cm~^ s~"^/LASS count in the 2 — 10 keV energy band appropriate for 
a power- law spectrum with photon index 1.7 for each fiux. We have re-calculated the fiuxes 
using our spectroscopy results for Sample 2 as well. 

The luminosities listed in Tables 3 & 4 are intrinsic luminosities, i.e., absorption- 
corrected luminosity using the fitted Nh. We show the absorption-corrected luminosities 
because we are more interested in the intrinsic AGN power in the later analysis. 
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3.2. AGNs with Special Considerations 

In this subsection, we summarize the spectral features of the AGNs for which the devia- 
tion from the simple power-law absorbed by neutral gas has a major impact on discussing the 
spectral statistics, intrinsic luminosities, and absorbing column densities. In particular, we 
have taken special attention to AGNs which gave unusually hard power- law index (F < 1.4) 
in the first step simple absorbed power-law fits, since they are likely to be heavily affected by 
partially covered and/or warm absorbers. Six AGNs in the entire sample have been found 
to fall in this category as detailed below. 

In calculating the intrinsic luminosity, we have used the spectral models involving partial 
and/or ionized absorbers. The global normahzation has been adjusted to give the observed 
HEAO-1 A2 (Al) count rate, if the AGN belongs to sample 1 or sample 2. The intrinsic 
luminosity in Tables 3 & 4 have been corrected for all absorbers. For the column density in 
these tables, we have used the sum of all neutral absorbers (including partial absorbers) as 
a representative value, but have not included those of any warm (ionized) absorber compo- 
nents. We explain our treatments for these 6 AGNs below. 

NGC 7582 A simple absorbed power-law fit to the ASCA data gave F = 1.2 and A^h = 
6 X 10^^ cm~^. The best description of its X-ray spectrum was given by Turner 
et al. (2000) from the Beppo-SAX data. Their model included an intrinsic power-law 
component of F = 1.95 absorbed by two components: a totally covering A^h,22 = 14.4 
absorber and a thicker absorber with A^h,22 = 160 blocking ~ 60% of the X-ray source. 
We employ this model for de-absorption (see above) and the sum of these two absorbing 
columns for the representative neutral column density. 

3A 0557-383 The hard X-ray source 3A 0557-383 was identified with a Seyfert 1 galaxy. 
Turner et al. (1996) reported the complex absorption below 2 keV of the X-ray spectrum 
from an ASCA observation. They found a good fit with either 96 % of the source 
covered by a column of 3.1 x 10^^ cm~^ of low-ionization gas or full covering by a 
column of 3.6 x 10^^ cm~^ of highly ionized gas, and an neutral iron K-shell emission 
line with equivalent width 300 eV. We used the latter model (Model D of Turner et al. 
1996). 

Kciz 102 The results of detailed analysis of our own ASCA data of this object have been 
reported by Miyaji et al. (2003) (M03). Kaz 102 has no sign of heavy absorption or 
a Fe Ka line. A single power-law with F ~ 1 described the ASCA spectrum well, 
while it had softer spectra when it was observed by ROSAT. This hard spectrum is 
probably due to extreme warm absorber component and/or refiection as discussed in 
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M03. In this study, we employ the warm absorber model with intrinsic F ~ 1.9 in M03 
for absorption correction. No neutral absorber is detected in the analysis and thus we 
assumed a zero A^h for the neutral gas. 

3C 445 The nearby broad-hne radio galaxy (BLRG) 3C 445 has very hard power law r=1.08 
and a large intrinsic absorption (A^h ~ 1-4 x 10^^ cm~^) in our first-step analysis of 
the XMM-Newton observation data. In Sambruna et al. (1998), the X-ray spectrum of 
3C 445 observed with ASCA were described by a dual absorber with an intrinsic flat 
power-law (F ~ 1.3) continuum. Wc have made our own analysis of the XMM-Newton 
data with two intrinsic absorbers, neutral partial covering and full covering cold gas in 
0.6 — 12 kcV band like Model C of Sambruna et al. (1998). Our analysis show a hard 
intrinsic power law (F =1.11±0.05) absorbed by partial covering (90 ±10%) neutral gas 
with A^H ~ 1.32±0.06 x lO^^ cm-^ and no full covering cold gas (x^/d.o.f = 1.27/159). 

We employ our XMM-Newton result for intrinsic luminosity calculation and statistics. 
Taking Sambruna et al. (1998) 's model had neghgible effects on the analysis in next 
sections. 

MKN 6 A detailed spectroscopic analysis of this object is given by Immlcr ct al. (2003), 
who analyzed the same XMM-Newton data combined with Beppo-SAX data. Their 
joint analysis preferred a model (Model 4) with a fixed-index power law (F=1.81) and 
two partially covering neutral absorbers (56±2 % of 8.1±0.5 x 10^^ cm~^ and 91±1 % 
of 2.3±0.1 X 10^^ cm~^), with a very small reflection component. We have re-analyzed 
the XMM-Newton using this model and checked the spectral fitting results, when F 
was set to free parameter. We found a slightly harder spectral index and small column 
density, but the difference of derived intrinsic X-ray luminosity is less than 10%. This 
has negligible effects on our HXLF and other statistical conclusions. We listed the sum 
of column densities of the dual absorbing gas and F = 1.81 and used the spectral shape 
of Immler et al. (2003) 's Model 4 for the absorption correction. 

NGC 4151 This Seyfert 1.5 galaxy NGC 4151 has long been known that it appears to have 
a flat power law (F ~ 1.3 — 1.7), a substantial column density of absorbing gas (Nh 
~ 10^^ cm~^), a strong soft excess at energies below 2 keV, and the broad iron K a 
emission line. We take the result of Model 2 in Schurch et al. (2003), who modeled 
the spectrum with two absorption components, a partially photoionized absorber and 
a neutral absorber, and a neutral reflection continuum. The column density of the 
neutral absorber is listed in Table 3 and the absorption correction has been made with 
their Model 2. 
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3.3. Other Spectral Features 

We have also examined the spectral fits in the broad band spectrum more in detail, 
in particular, paying attention to the existence of iron lines and/or a soft excess. We have 
started with a power-law with the photon index fixed to that in first step analysis, with the 
Galactic absorption at the AGN position. For those where the simple absorbed power-law 
model have significant residuals in the broad-band spectrum, we considered more complicated 
models, involving, e.g., a soft excess and or an iron line. The fitting results with the soft 
excess are shown in Table 5. We listed the only best fit results for soft excess, within several 
spectral model fittings (blackbody, broken-power-law, 2-powerlaw, soft thermal emission, 
dual absorber, etc.. ). Typical spectrum of the sample in our spectral fitting are shown 
in Figure 3. For those the addition of an iron line improved the fit, the central energy, 
normalization, and equivalent width the line, represented by a single Gaussian, are listed in 
Table 6. In the spectrum of MKN464, there are significant residuals in 1.5 — 3.0 keV band 
and this range has been excluded for the analysis. 

The photon index and intrinsic column density for most of all AGNs are consistent to 
that in first step analysis within errors, even if the photon index and the intrinsic absorption 
(A^H in "zwabs") are set to free parameters in this second-step spectral analysis. Thus the 
existence of soft excess and iron emission line does not have significant effect in our further 
analysis. 

The statistics of such features has its own importance. However, in our current data, 
the limits on the detection of these features vary observation by observation and difficult to 
model. Thus one should not take advantage of the good completeness of our sample to make 
a statements on the statistical properties of the iron lines and soft excess in combination 
with our sample definition. 

4. STATISTICS OF SPECTRAL PROPERTIES 

4.1. Neutral A^h and Photon Index Distribution 

Figure 4 shows the distribution of the absorbing column densities (Nh) of neutral ab- 
sorbing gas intrinsic to the AGNs. As described in Sect. 3.2, we have added the value of 
any partial covering neutral column densities to the Ah value used here. On the other hand, 
the column densities of ionized absorbing gas have not been included (3A 0557, Kaz 102, 
NGC 4151). This is because the neutral gas and ionized gas come from different regions, the 
former being in a molecular torus at about 10^^ — 10^° cm from the nucleus, while the latter 
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at a region closer to the nucleus (~ 10^^ cm, (e.g. Reynolds & Fabian (1995)). Thus they 
do not have much relevance to each other. 

We also plot an expected relative spatial number density in each A^h bin, assuming a 
Eucleadian log N — log S relation for the intrinsic flux distribution of the X-ray AGNs: 



in which, A^spatiai is spatial number density, /intrinsic is intrinsic flux of primary X-ray com- 
ponent, /absorbed is observcd flux (after absorption), and A^raw is the raw number in the A^h 
histogram. Upon calculating the /intrinsic of each bin, we have assumed that the intrinsic 
power-law source of F = 1.7 was absorbed by cold gas of mean column density of each bin. 
This may not be correct in case of the source has partially covering absorbers. However, the 
effect is much smaller than the errors due to small number statistics in each bin. 

Figure 5 shows the distribution of spectral photon index. The F values of NGC 3227, 
3A 0557, Kaz 102, H 1537+339, and MKN 6 are excluded from the histogram, because they 
have been flxed during the spectral flt, or no detailed spectral information. The mean photon 
index is 1.76 with a — 0.2. This is similar to the result of WiUiams et al. (1992), who studied 
an incomplete sample of 13 QSOs observed by GINGA and obtained the result of the mean 
photon index of 1.81±0.31. Even if we use the results from flrst-step single absorbed power- 
law flt results (Sect. 3) the 6 sources discussed in Sect. 3.2, the mean photon index has not 
change remarkably, giving 1.72±0.24. In George et al. (2000), the mean photon index has 
~ 2.0 for various spectral analysis, while the relatively wider distribution (1.3 < F < 3.0). 
And in Reeves & Turner (2000), they derived mean photon index 1.89 ± 0.05 with a = 0.27. 
The two other observations are somewhat softer distributions, which is assumed be due to 
the effect of a reflection component. We have found no correlation between F and hard X-ray 
luminosity, which is consistent with results of George et al. (2000), and Reeves & Turner 



4.2. Relation between Hcird X-ray Intrinsic Luminosity and A^h 

Figure 6 shows the scatter of the AGNs, in the hard X-ray luminosity - intrinsic neutral 
A^H plane. Different symbols are used for AGNs with different optical types, i.e. type 1 
(Seyfert 1-1.5, typel QSO, or BLRG), type 2 (Seyfert types above 1.8) AGNs or LINER 
candidate. Most of optical type 1 AGNs have small X-ray absorption, log A^h < 21.5, and 
most of Seyfert 2 galaxies have large X-ray absorption, log A^h ^ 21.5 as widely known. 
Some AGNs optically classifled as type 1 have large A^h (log A^h ^ 21.5) values. They 




(1) 



(2000). 
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usually have X-ray spectra involving partially covering absorbers. We have include such 
sources in the absorbed AGN group for the HXLF analysis. The optical classifications in 
this paper are not necessarily accurate, because they are based on the detection significance of 
broad emission lines and depend on the quality of the optical spectra. Indeed, there are some 
AGNs (ESO103-G35, NGC526a, 3c445, etc..) for which different sources give different optical 
classifications. However, our spectral results may suggest that the anomalous cases reflect 
the existence of substantial variance in the dust to gas ratio, or a geometrical separation to 
the line of sight. More significant data analysis with careful regard to the variability would 
derived whether these objects are common in AGNs or not. 

We observe no AGN with log > 44.5 and log Ah ^ 21.5. This is consistent with 
U03's analysis including ASCA and Chandra surveys, which find a decrease of absorbed AGN 
fraction towards high intrinsic luminosities. A further discussion of this effect will be given 
in Sect. 6. 



5. LOCAL HARD X-RAY LUMINOSITY FUNCTIONS 

Our purpose of this section is to calculate the local Hard (2 — 10 keV) X-ray Luminos- 
ity Function (HXLF) of AGNs from our sample. We have taken the U03's approach and 
constructed the 2 — 10 keV HXLFs for the intrinsic X-ray luminosity. We have divided 
the sample into unabsorbed and absorbed AGNs at log A^h = 21.5 for the use with HXLF 
analysis. The number of X-ray absorbed AGNs is 16, while that of X-ray unabsorbed AGNs 
is 33 in our sample. 

We have used the ^V~^ estimator (Avni & Bahcall 1980), which is a modified version 
of original ^V^^ (Schmidt 1968). We can write the binned HXLF as, 

-(logL.)^S. ^y^% , (2) 



dlogL/ (AlogL.) 

where is the available comoving volume in which the i-th AGN would be in the sample, 
within the redshift range. The value can be expressed by 

dA{zfc—{z) ■ A,s{CR) dz (3) 

-min 

where dxiz) is the angular distance, dT/dz{z) is the differential look back time per unit z, 
CR is the observed count rate of the AGN if it was placed at the redshift z, and ^eff(CR) is 
the survey area as a function of the count rate. We take ^min = and ^max = 0.4 to include 
all the AGNs in the sample. 
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Since our sample consists of two subsamples defined in the count rate limits from two 
different detectors, our 1^, consists of two parts: 



where T^*^ and comoving volumes of the universe within which the i-th AGN 

would be observable in the Sample 1 and Sample 2 respectively. In calculating V^2 '^^ AGN 
in Sample 1, we have used the predicted LASS cts cm~^ of HEAO-1 Al experiment, 
which has been estimated by using the HEAO-1 A2 R15CR, the spectral shape of this 
source, and the detector responses of HEAO-1 A2 experiment, and vice versa. In the HXLF 
calculation, we have used the observed flux (i.e. after absorption) in estimating the V^, while 
X-ray intrinsic luminosity (de-absorbed luminosity) has been used to divide the sources into 
luminosity bins. 

As mentioned in Sect. 2, we have used the effective survey area for Sample 2 after 
correcting the LASS count rate for the effects of confusion noise. This correction is discussed 
in detail in Appendix A. After the correction, the area as a function of the corrected count 
rate gradually decreases as the CRcon decreases to the hmit of 1.9 xlO~^ LASS cts cm~^. 

A small number of very luminous AGNs in our sample are observable at redshifts where 
evolution has a non-negligible effect {z > 0.2 — 0.3). Thus we have also derived the evolution- 
corrected HXLF by giving weights to the calculation. 



with which we can estimated the HXLF at z = Zq by including the density evolution factor 
[(1 + 2;)/(l + 2;o)]*'. We use the index p = 5, which is consistent with the soft X-ray luminosity 
function evolution by Hasinger et al. (2005). The impact of the evolution to the XLF was 
found to be as small as ~ 10% in the 45.5 < log < 46.5 bin. Since our sample is mainly 
local, K-correction has a negligible effect. Most luminous AGNs are unasborbed with the 
photon index ~ 1.9, with very little K-correction even at z — 0.4. Absorbed AGNs typically 
have such low luminosities that they are not detectable at high enough redshifts for the K- 
correction to be important. The K-corrections are typically < 5% for the absorbed AGNs. 

Figure 7 shows the Hard X-ray luminosity function calculated from our sample for the 
evolution-corrected case. We take 2^0 = for Figure 7. We have made the luminosity 
bins same in all HXLF calculations for comparisons. We only have upper limits for absorbed 
AGNs in the three highest luminosity bins and we plot only interesting (90%) upper limits, i.e. 
those which make significant constraints. Figure 7 clearly shows different X-ray luminosity 
function slopes between absorbed and uabsorbed AGNs at the high luminosity end. The 
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ratio of absorbed / unabsorbed XLF is 2.35±2.1 in log Lx ~ 42 - 44, and 0.11±0.08 in 
44 — 44.5). The significance of the difference of the absorbed/unabsorbed ratio in these two 
luminosity bins has been estimated using a bootstrap resampling. Out of all the 41 AGNs 
in the log Lx ~ 42 — 44.5 range, we have randomly resampled the same number of AGNs 
allowing duplicates. The ratio of the absorbed / unabsorbed HXLFs in the two luminosity 
ranges have been calculated for each of 10000 bootstrap-resampled samples. As a result, the 
probability that the ratio reverses, i.e., absorbed / unabsorbed ratio becomes greater in the 
higher luminosity bin, turns out to be only less than 0.4%. 

In order to compare absorbed and unabsorbed AGN HXLFs more quantitatively, we also 
made a maximum likelihood (ML) technique. The ML-fitting of the HXLF to the data where 
individual AGNs have different spectra is not straightforward. We use U03's approach in 
fitting A^H function and the intrinsic HXLF. There is a limitation in this approach, however, 
that we have to assume a simple absorbed power-law spectrum for each source. For the 
fits, we have used the sources in the range z < 0.2 and log A"h < 24 to compare with the 
U03 results. We again divided into unabsorbed and absorbed AGNs at log Ah = 21.5 in 
this sample. In this analysis, the Ah function was assumed to be constant within each of 
log A^H < 21.5 and 21.5 < log Ah < 24 bins. The fits do not depend significantly on details 
of the evolution within the small redshift range. Thus, a pure density evolution model by 
U03 {e{z) = (1 + zy-"^) was used to represent the small amount of evolution. The best-fit 
parameters of the XLFs were calculated following the same procedure as U03. 

Table 7 shows the best fit parameters from the ML-fits. The fit result from U03's PDE 
model is also listed. The two power-law indices for the 'AH' sample are consistent with that 
in U03, while the normahzation A seems larger than the U03 value. This is partially because 
A is strongly depend on the log L*, parameter, for which we found a somewhat smaller 
value. Also, U03 is based on the global fit to all redshifts, where the higher redshift samples 
contribute much of the constraints to the relevant parameters. In comparing absorbed and 
unabsorbed AGN samples, the ML fits have been made with a fixed log = 44.0 in order 
to avoid strong dependence between fitting parameters. This value is consistent with the 
best-fit values for both samples, when fitted with this parameter free. As expected, the 72 
of absorbed sample is significantly steeper than that of unabsorbed AGNs. 

In Fig. 7 (evolution-corrected to z = 0.1), we overplot the best-fit result of smoothed two 
power law formula for SXLF in the < 2; < 0.2 bin by Hasinger et al. (2005). The 0.5 — 2 keV 
SXLF has been converted to 2 — 10 kcV band using a F = 1.8 power-law spectrum. Our 
evolution corrected HXLF of unabsorbed AGNs obtained from our sample is marginally 
consistent with the Hasinger et al. (2005) best-fit model. We also overplot the results from 
Sazonov (2004) for unasborbed and total AGNs. They used the 3 — 20 keV band, where the 
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attenuation of flux by absorption have only minor effect (by ~ 8% for log A^h =22.5). We 
have thus converted their 3 — 20 keV HXLFs (for luminosities uncorrected for absorption) 
to 2 — 10 keV using a F = 1.8 power-law for overplots. Detailed comparisons are discussed 
in Sect. 6. 

6. DISCUSSION 

6.1. Comparisons of Luminosity Functions and Volume Emissivity 

We have calculated the local HXLF with binned method for our sample defined 

with HEAO-1 experiment. The HXLF for our log N-^ < 21.5 results is marginally consistent 
with the best-fit model of SXLF in the < 2; < 0.2 bin in Hasinger et al. (2005) as well as with 
the results from Sazonov (2004) for unasborbed and total AGNs. It is not surprising that 
we have small discrepancies because our log A^h < 21.5 selection does not necessarily match 
with Hasinger's soft X-ray type 1 selection. A notable exception is the 44 < log Lx < 44.5 
bin, where our sample gives an excess which cannot be explained by statistics alone. 

The local X-ray volume emissivity of brightest AGNs can be estimated by Ej(Lx,i/l2). 
The volume emissivity value estimated from our sample is P2-10 keV — 5.85 ± 1.17 x 10^^ erg 
s"-*^ Mpc""^, which is derived using observed luminosity (absorbed). Note that there is no 
AGN with log Lx < 42 erg s"^ in our sample. 

Miyaji et al. (1994) estimated the 2 — 10 keV local volume emissivity by modeling the 
cross-correlation between the all-sky cosmic X-ray background (CXRB) surface brightness 
from the HEAO-1 A-2 experiment and the IRAS galaxies, obtaining (6.1 ± 1.7) x 10^^ erg 
s~^ Mpc~^ (for Hq — 70 km s~^ Mpc^^). This is consistent with the AGN volume emissivity 
from our sample, implying that the contribution of lower luminosity sources (e.g., star- 
forming galaxies and LINERS) as well as other sources, non-active galaxies to the local 
volume emissivity is very small. On the other hand, Sazonov (2004) estimated the X-ray 
(3 — 20 keV) volume emissivity of log Lx > 41 erg s"-*^ from emission line AGNs with RXTE 
All Sky Slew Survey and found p3_2o kev(> 41) = (4.9 ±1.2) x 10^^ erg s"^ Mpc-^ (for Hq = 
70 km s~^ Mpc~^ z > 0.001), which can be converted to (2.7 ± 0.7) x 10^^ erg s~^ Mpc~^ 
in the 2 — 10 keV band using the absorbed power-law model of AGNs estimated based on 
the RXTE hardness ratios. Their local volume emissivity estimate is only about a half of 
our current estimate. One possibility is an effect of local large-scale structure z < 0.015, 
where our sample and the RXTE sample have different weights. 

The HXLF from our sample seems larger than that of RXTE with log L^r^ AA — 44.5. 
There is a possibility of some systematic in Sazonov (2004) 's 2 — 30 keV to 2 — 10 keV 
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conversion. Their absorbed F = 1.8 power-law model fit to the RXTE hardness ratio based 
on 3 — 8 and 8 — 20 keV count rates is only sensitive to log A^h > 22. Also their harder band 
is likely to be affected by reflection bumps. At this moment, we are not certain whether the 
differences in the volume cmissivity and HXLF results are due to the large-scale structure 
(cosmic variance) effect, assumed spectra in band conversions, or other systematics. A 
similar study with a large-scale sensitive hard X-ray survey with sufficient AGNs to exclude 
the z < 0.015 universe and to sample a large volume of the present-day universe, such as 
available from eROSITA or MAXI will give an ultimate solution. 



6.2. Luminosity — Absorbed AGN fraction Relation 

One way of testing the simplest unified scheme, where the difference between absorbed 
and unabsorbed AGNs are purely from the orientation effect and the intrinsic properties are 
essentially the same is to compare the shapes of unabsorbed and absorbed AGN XLFs in 
intrinsic (i.e. de-absorbed) luminosity. We have found a difference between the HXLFs of the 
X-ray absorbed AGNs and that of unabsorbed AGNs in that the number density of absorbed 
AGNs drops more rapidly than unasborbed ones towards high luminosities. This trend can 
be also seen in the scatter of the sample AGNs in the hard X-ray intrinsic luminosity versus 
intrinsic neutral A^h diagram in Figure 6, in which we see a void of AGN in log Lx > 44.5, 
log A^H > 21.5. This trend has been clearly demonstrated by U03 as the decreasing absorbed 
fraction of AGNs as a function of intrinsic luminosity from a larger sample including AGNs in 
this work, from ASCA and from Chandra surveys. This lead U03 to suggest a modified AGN 
unified scheme, where the difference between type 1 (unabsorbed) and type 2 (absorbed) 
AGNs is not solely due to the viewing angle effect, but involves some intrinsic difference of 
the geometric and/or physical conditions of the absorbing tori around the SMBH based on 
the radiative power of the primary X-ray component. 

To verify that the high-luminosity high-absorption void really reflects the actual de- 
ficiency of high-luminosity highly absorbed AGNs rather than by a selection effect in our 
sample, we calculate the expected number of AGNs in the case where the simplest unified 
scheme was true, i.e., there was no intrinsic difference in the absorbing geometry as a func- 
tion of AGN power. We have started with the smoothed two-power-law expression of the 
0.5 — 2 keV soft X-ray luminosity function in the 0.015 < z < 0.2 bin by Hasinger et al. 
(2005). We have slightly tweaked the parameters and added AGNs in the absorbing columns 
of log A^H=20, 21., 22, 23, and 24 with the same luminosity function shape, allowing the 
relative normalization to vary, following an approach similar to an AGN population synthe- 
sis model. The ratios of number densities in each A^h bin have been adjusted to match the 
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actual numbers of AGNs in log A^h ranges of < 20.5, 20.5 - 21.5, 21.5 - 22.5, 22.5 - 23.5, 
23.5 — 24.5 respectively. A further tweak of the parameters have also made to reproduce 
the number of AGNs in log Ah < 21.5, log Lx > 44.5 region (luminous unasborbed AGNs). 
An important assumption in this model is that the shape (all parameters besides normal- 
ization) are identical for all A^h classes. The expected number of AGNs in the "luminous 
absorbed" void under this model turns out to be five, while we observe none. Four of the 
five AGNs expected under this model would be in Sample 1, where only one of 57 X-ray 
sources are unidentified. Treister et al. (2005) suggested that apparent such effect could 
arise from selection effects, their analysis is based on a Chandra survey focusing on a much 
fainter population with limited identification completeness. Their argument does not apply 
to our analysis or that of U03, involving highly complete samples. (See also Treister & Urry 
(2005), where they also found this effect by including more samples.) 

Recently Zhang (2005) made an interesting suggestion that, if most of the intrinsic 
power-law component in the X-ray spectra of AGNs is from a Comptonized emission by 
optically thick hot plasma above the accretion disk, this relation can be explained solely by 
the viewing angle effect due to geometrical projection. This assumption needs a physical 
condition that the planes of the accretion disk and the dusty torus are co-aligned. We note 
also that the Comptonizing plasma need to have a planer geometry and to be parallel to the 
accretion disk. If the emission region is optically thin and/or composed of hot spots rather 
than a planer region, the projection effect would not arise. 



6.3. Flux Variability 

In our investigation, a complete flux-limited sample of brightest objects were re-observed 
with XMM-Newton, ASCA and/or Beppo-SAX. We plotted the re-observed flux versus 
HEAO-1 observed flux in Fig. 8 for unabsorbed and absorbed AGNs. The HEAO-1 flux 
has been corrected for the confusion noise bias described in Appendix A. We see from Fig. 8 
that most AGNs show lower fluxes upon the re-observations, even after the correction. 

The decrease can be naturally explained as a result of AGN variability and the log — 
log S distribution of the AGNs which rises rapidly as flux becomes lower. A flux-hmited 
brightest sample is more hkely to pick up an AGN's brighter than average state, while the 
re-observation picks up all states equally. We formulate this effect as follows. 

At sufficiently bright fluxes, the log N — log S relation follows the Euclidean relation 



N{> S) oc S'^-^, 



(6) 
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or in the differential form 

n(S) = dN(S) /dS (X 5-2-5 (7) 
For the convenience of the calculations, we define I = ln S. Then, Eq. 7 can be rewritten as: 

n{l)dl oc e-^-'^^dl (8) 

Suppose that the variability of each individual AGN is characterized by the Gaussian 
log-flux distribution with a standard deviation a: 

v{l',l) = (27ra2)-V2exp[-(/ - I'Y/i^a^ (9) 

the apparent source count above a log flux of Iq at the time of the survey observation can be 
expressed by 

/oo poo 
{n{l') / v{l',l)dl}dl', (10) 
oo J lo 

If an AGN whose average log flux is I' was observed at I at the time of the survey, the 
expectation value of its flux upon the re-observation is I'. Therefore the expectation value 
of the decrease of the flux {{dl) = {I' — I)) can be expressed by 

...._ mninj,:vii':i){i-ndi}di' 

IZ.{n{ni:v{l':l)dl}dU 

Under the assumption that n{l') follows the Eucleadean relation Eq. 8 and log normal 
flux distribution of AGN flux variation in Eq. 9, and by further assuming that all AGNs 
vary at a typical amplitude characterized by a, Eq. 11 can be integrated numerically. Figure 
9 shows the numerical solution of Eq. 11 in the a — (dl) plane. /,From Fig. 9, we see that (/) 
is always positive, which means re-observation give lower fluxes on average. This is partially 
because there are much more fainter AGNs under the distribution Eq. 7. Thus the net effect 
of AGN variability is that AGN in the flux-limited sample seems to show a flux decrease 
upon the re-observation on average. In turn, from the average decrease of the flux, we can 
estimate the typical AGN variability amplitude. In our sample, the re-observed flux by 
ASCA, XMM-Newton or Beppo-SAX was smaller than the HEAO-1 flux by a factor of 0.28 
on average (after correction for spectral response and confusion). Figure 9 shows that this 
corresponds to cr = 0.914 or a typical variability amplitude of a factor of 2.5. 



7. CONCLUSIONS AND PROSPECTS 



We summarize our investigation and main conclusions below. 
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1. We have defined a complete flux- limited sample from a combination of HEAO-1 A2 and 
A1/A3 surveys with a well-defined flux-survey area relation. The sample is composed 
of 28 AGNs from Piccinotti et al. (1982) and 21 AGNs from a deeper sample by Grossan 
(1992). 

2. A bias to the measured flux due to confusion noise have been found to be signiflcant 
in Sample 2 (A1/A3). We have modeled and made flrst-order corrections to the fluxes 
and survey areas. 

3. We have obtained X-ray spectral information for all (but one) AGNs in the combined 
sample from XMM-Newton, ASCA and Beppo-SAX. 

4. The spectra have been first modeled as an absorbed power-law. We have made closer 
look at 6 AGNs which gave spectral index of F < 1.4 in this fit and used models 
involving (multiple) partial covering and/or ionized absorbers. 

5. We have found the mean photon spectral index of F 1.76 with a la dispersion of 
0.2. 

6. We find the A^h distribution for our sample, which can be used as a constraint to X-ray 
population synthesis modeling. 

7. We have constructed local hard X-ray luminosity functions (HXLFs) from our sample, 
separately for absorbed and unabsorbed AGNs, as well as for both, as a function of the 
intrinsic 2 — 10 keV luminosity. We have also made HXLFs which have been corrected 
for density evolution. 

8. The scatter of the sample AGNs in the versus intrinsic luminosity plane shows 
a void of AGNs at a high-luminosity, high absorption regime, where we expect ~ 5 
AGNs if the intrinsic luminosity functions between absorbed and unabsorbed AGNs 
were identical. This, as well as the difference in the HXLFs between absorbed and 
unabsorbed AGNs, is not likely due to selection effects. 

9. The X-ray fiuxes of AGNs in the sample observed by XMM-Newton, ASCA and/or 
Beppo-SAX were on average lower than those observed by HEAO-1 for the same ob- 
jects. The mean flux decrease upon re-observation of AGNs in a flux-limited sample 
has been formulated can be naturally explained by the AGN variability, where the 
AGNs are hkely to have brighter-than- average flux at the time of the sample-deflning 
survey. 
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As a concluding remark, we note that the sensitivity of large-area surveys at hard X-ray 
bands currently available is still not sufficient, in contrast with the soft band, where ROSAT 
All-Sky Survey produced enormous dataset to sample the present-day universe unabsorbed 
AGN population. In this work, we had to depend on data obtained in the 1970's to define 
our sample, to include the absorbed AGN population. Because of this, our investigation is 
subject to fundamental limitations in terms of object count statistics, confusion and mis- 
identification problems. To obtain sufficient number of AGNs, we also had to use very nearby 
objects {z < 0.015), where the quantities such as HXLF or volume emissivity is subject to 
the density fiuctuations due to the local large scale structure. The understanding of AGN 
evolution is never complete until we sample the present-day universe fairly. In the near 
future, a slew survey from the Swift BAT mission (Markwardt et al. 2005) and a dedicated 
all-sky exposure with INTEGRAL will produced a less biased bright end catalog of AGNs and 
after a spectroscopic followup such as this work, we expect improvements in our knowledge 
in this region. Future missions such as MAXI^ and eROSITA would find numerous AGNs in 
our local universe will enable us to make a precision statistics of AGNs in the present-day 
universe. 
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A. Correcting for Confusion Effects on HEAO-1 Al Counts 

Our sample selections are partially based on count rates (CRs) observed by the HEAO- 
1 Al experiment, which have been cataloged by Wood et al. (1984) (hereafter W84). The 
CR limit imposed by Grossan (1992) of C-Rum = 3.6 x 10"^ LASS cts s^^ cm^^ is subject 
to confusion noise, in addition to the formal error given by W84, which only includes the 
effects of photon counting statistics. Due to these effects, the CRs are subject to systematic 
overestimations near the detection limit. This may bias the estimates of luminosity function 
derived in Sect. 7. 

We have made a first-order correction to these effects by Monte-Carlo simulations. We 
have taken the following steps. 

1. Sources are generated based on approximate log — log S relation, extrapolated down 
to approximately one order of magnitude below our flux limit. In this simulation, 
LASS CRs (Ci?Ai) have been used rather than fluxes in physical units, because we 
define our limit by CRs. From the W84 catalog, we have estimated: 

N{> CRai) = 7 X 10-\CRai/5 x 10-2)-^-^[deg-2] (Al) 

2. For each generated source with true underlying CR (Ci?gen), an observed CR (CRohs) 
has been simulated by 

CRohs = CRgcn + ACi?err + ^C R^s, (A2) 

where ACi?crr simulates the deviation due to photon counting statistics and AC Rets 
due to confusion noise as detailed below. 

• For a given CRgen, an error a is selected from objects in |6| > 30, \EB\ > 35, 
and CR < 0.01, where EB is the ecliptic latitude. ACi?err is derived from a 
random Gaussian deviation using the selected a. This may underestimate errors 
for brighter sources, but the effects we consider is only non-negligible for sources 
near detection limits. 

• ACRcfs is derived from Monte-Carlo simulations, where objects are generated 
based on Eq. Al in the range 5 x 10~^ < CR < 3.6 x 10^^ (the upper bound is 
our flux hmit) and distributed over the HEAO-1 Al beam (W84). The total CRs 
folded with the beam (Ci?beam) have been calculated for 2000 runs. If C Rheam is 
greater than 3.6x10"^ above the mean value, the run is rejected. 

• The mean value ((Ci?beam)) has been recalculated after the rejection. Then 
ACRcis has been calculated as Ci?beam — (C'-Rbeam), where CRheam is taken from 
one of the simulation runs which was not rejected. 
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3. In order to simulate detected sources, only those with CRohs > 3.6 x 10~^ have been 
selected (C-Robs of the selected objects is indicated by CR^ct, representing 'detected 
count rates'). Figure 10(a) shows the scatter diagram of CR^ct versus CRgen for those 
objects. The simulated objects have been then binned by CR^ct and the mean value 
of CRgen {{CRgcn}) has bccu calculated in each bin. This gives a measure of mean 
expected true CR of the sources as a function of the detected CR. 

4. Our next step is to find a smooth analytical representation of {C Rgcu) / C R^et as a 
function of CRobs- Since we know that {C Rgcn) / C R^et goes asymptotically close to 
unity at higher count rates, we have used the form: 

Ccr = (Ci?det)/C-Rdet = 1 " exp[-(Ci?det " Scr)/«'cr)], 

{CRdet > CR) (A3) 

where Ccr represents the CR correction factor. By a fit with the standard deviation 
of CRgcn in individual bins as a, we find the best-fit values: 

Scr = 1.6 X 10"^; Wcr = 2.6 X 10"^ (A4) 

The binned {C Rgcn) / C Rdct values are plotted against CRdet in Fig. 10(b) along with 
the best-fit function form. Note that we use with these cr's only as a weighting 
scheme for the fit and the formal parameter errors of the statistics does not have 
any meaning. 

5. The detected CRs have been corrected by CRcow = e'er CR^^i. By this correction, 
the estimated 'true' count rate for a source detected at the faintest limit of CR^et — 
CRnn, = 3.6 X 10"^ becomes Ci?corr,iim = 1.9 X 10"^ 

6. In calculating the XLF, we also have to correct the survey area for incompleteness 
of the detection, where not all sources with a true count rate of, e.g. > 1.9 x 10~^ 
will have CRobs > 3.6 x 10~^. This is made by dividing the differential logA^ — logS* 
curves of CRcon by that of CRgcn calculated using the simple geometrical area (Ageo) 
of the survey. Using the same functional form as Eq. A3, the effective survey area as 
a function of Ci?corr is expressed by 

^eflfMgeo — ^ — exp[— (Cit!corr " Sa)/wa], (C'it^corr > C'il!corr,lim) 

— 0, {CRcoii < C'i?corr,lim) ■ i-^^) 

We have found the best-fit values of the parameters: 

SA = -3.3 X 10"^; WA = 8.5 x 10"^ (A6) 
The resulting effective area curve is plotted against Ci?corr in Fig 10(c). 
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7. Using the corrected Ci?corr and Aes, a corrected log A^ — log S curve has been calculated 
to see that it is in reasonable agreement with the initially assumed one. This curve 
is plotted in Fig 10(c) along with those from CR^et and Ci?corr using the geometrical 
area. 

8. Now we are ready to calculate the luminosity function from the real data by applying 
the corrections Eqs. A3 & A5. to the cataloged CR and geometrical area respectively. 

We have also made the same experiment for Sample 1 based on the HEAO-1 A2 Pic- 
cinotti et al. (1982) catalog. The effects are negligible in the case of Sample 1. 
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Galactic Longitude 



Fig. 1. — Spatial distribution of our sample. Dashed line represents b=±20 degree, dotted 
line represents the selected region of Sample 2 (Grossan). Open circle : Sample l(Piccinotti), 
Cross : Sample 2 (Grossan). 



-29- 



47 



46 



X X 



I 

o 



45 



bo y 



.X 



X 



> 

o 



44 



CCD 



X. 



X_. 



nX 



M 



X 



X 



X 

o 

h-1 



43 



42 



41 



10 



10 



10 



Redshift 



Fig. 2. — Relation between X-ray luminosity and redshift of the sample. Open circle : Sam- 
ple 1 (Piccinotti). Cross : Sample 2 (Grossan). For Sample 2, the flux is calculated from the 
corrected count rate described in Appendix A. 
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Fig. 3. — Examples of spectral fit results in the E I{E) scale, for IIIZw2, NGC 7314, 
ESO 103-G35, and Fairall 9 in 0.3-10 keV band. The spectra are represented by thin red 
crosses, and green and blue lines represent the all model spctrum and each spectral compo- 
nent, respectively. The details of each spectral component are listed in Table 3-6. Figures 
3.1 - 3.5 are available in the electronic edition of the Journal. 
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Fig. 4. — The intrinsic Nh (neutral) distribution of the AGNs in our sample. We also plotted 
spatial number density calculated from Eq. 1. 
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Fig. 5. — Photon index distribution of the sample. The open histogram represents all the 
sample, the filled histogram represents X-ray absorbed AGNs. 
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Fig. 6. — Relation between hard X-ray luminosity and intrinsic Nh- Type 1 galaxy (Seyfert 
1, 1.5, Type 1 QSOs) are shown as filled circles and Type 2 AGNs (Seyfert 1.9 - 2) are 
shown as open circles respectively. One probable LINER is plotted with a filled star. The 
luminosity is intrinsic (before absorption). 
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Fig. 7. — Hard X-ray luminosity functions for all AGNs (black, square), X-ray unabsorbed 
AGNs (red, circle), and absorbed AGNs (blue, triangle) in our sample for the evolution- 
corrected case. Thick solid lines and dotted lines show the ML-fitting results of the XLFs for 
unabsorbed and absorbed AGNs, respectively. Thin solid lines show type 1 AGN XLF from 
Hasinger et al. (2005). Thin dashed and dotted lines show the results from Sazonov (2004) 
for all AGNs and unabsorbed AGNs respectively. All overplotted lines have been converted 
to 2-10 keV from the originals assuming a F = 1.8 power-law. Errors are la. 




Fig. 8. — Relation between ASCA or XMM-Newton observation flux and HEAO-1 observa- 
tion flux. Circle : unabsorbed AGNs. Triangle : absorbed AGNs. 
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Fig. 9. — The numerical solution of Eq. 11 in the o — {dl) plane, where a is typical deviation 
of the variability of AGN flux, and {dl) = (/' — /) represents the expectation value of the 
decrease of the flux. We can estimate a typical a of our sample from this a - (dl) relation 
under the assumption discussed in Sect. 6.3. 
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Fig. 10. — (a) The scatter diagram of the generated "true" count rate (CRgcn) versus the 
detected count rate (C-Rdct) in our simulation. The sohd hne shows CRgcn = CRdct- (b) 
The {C Rgen) / C Rdet value is plotted as a function of CR^et- The error bars indicate the la 
dispersion for the simulated object in the CR^et bin and not the standard deviation of the 
mean. The solid curve shows the fit to Eq. A3, (c) The corrected area curve for Sample 2 
from Eqs. A5 & A6. (d) The logA^ — logS curves for the simulations. The dashed hne 
shows the input model. Sohd triangles show the logA^ — log 5" calculated from CR^^t with 
the geometrical area. Sold squares show the case Ci?corr is used with the geometrical area. 
Open squares show the case after correcting both CR and the area (Ci^corr used with A^q). 
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Table 1. Source List - Sample 1 



Source ID 


Name 


RA 


Dec 




Redshift 


Galactic Njj 






(J2000.0) 


(J2000.0) 






(lO^O cm-2) 


1 


IC 4329A 


13 49 19.2 


-30 18 34 


SI 


0.0168 


4.42 


2 


3C 273 


12 29 06.7 


02 03 08 


QSO 


0.158 


1.79 


3 


NGC 2992 


09 45 41.9 


-14 19 15 


S1.9 


0.0076 


5.37 


4 


NGC 5506 


14 13 14.9 


-03 12 27 


SI. 9 


0.0072 


3.81 


5 


NGC 526A 


01 23 54.8 


-35 03 55 


SI. 5 


0.0202 


2.20 


6 


NGC 7582 


23 18 23.0 


-42 22 11 


S2 


0.0044 


1.59 


7 


ESQ 198-G024 


02 38 19.6 


-52 11 34 


SI 


0.045 


3.13 


8 


MCG -6-30-15 


13 35 53.3 


-34 17 48 


SI 


0.0082 


4.08 


9 


MKN 509 


20 44 09.0 


-10 43 15 


SI 


0.0352 


4.11 


10 


3C 120 


04 33 11.0 


05 21 15 


SI 


0.033 


11.1 


11 


NGC 7172 


22 02 01.7 


-31 52 18 


S2 


0.0085 


1.65 


12 


NGC 3783 


11 39 02.0 


-37 44 19 


SI 


0.0096 


8.50 


13 


MKN 926 


23 04 43.4 


-08 41 08 


S1.5 


0.047 


3.60 


14 


NGC 7469 


23 03 16.0 


08 52 26 


SI 


0.0161 


4.87 


15 


NGC 4593 


12 39 39.0 


-05 20 39 


SI 


0.0099 


2.31 


16 


NGC 3227 


10 23 30.5 


19 51 55 


S1.5 


0.0052 


2.36 


17 


ESQ 141-G55 


19 21 14.2 


-58 40 12 


SI 


0.0371 


5.09 


18 


IIIZw2 


00 10 31.0 


10 58 30 


SI 


0.0898 


5.72 


19 


Fairall 49 


18 36 58.4 


-59 24 07 


S2 


0.020 


7.13 


20 


3C 445 


22 23 49.0 


-02 06 12 


SI 


0.0562 


5.01 


21 


MKN 1152 


01 13 50.1 


-14 50 44 


S1.5 


0.0527 


1.69 


22 


NGC 7314 


22 35 46.0 


-26 03 02 


S1.9 


0.0046 


1.46 


23 


MCG -1-24-12 


09 20 46.3 


-08 03 22 


S2 


0.0198 


3.56 


24 


ESQ 103-G35 


18 38 20.3 


-65 25 42 


S1.9 


0.0133 


7.64 


25 


3A 0557-383 


05 58 02.1 


-38 20 05 


SI 


0.0344 


3.98 


26 


MRK 590 


02 14 33.6 


00 46 00 


SI 


0.0263 


2.68 


27 


H 1846-786 


18 47 02.8 


-78 31 50 


SI 


0.076 


9.05 


28 


Fairall 9 


01 23 46.0 


-58 48 21 


SI 


0.046 


3.19 



''SI : Seyfert 1, S2 : Seyfert 2 , S1.5 : Seyfert 1.5, S1.9 : Seyfert 1.9 
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Table 2. Source List - Sample 2 



Source ID 


Name 


RA 


Dec 


Type"* 


Redshift 


Galactic Nh 






(J2000.0) 


(J2000.0) 




(^) 


(1020 cm-2) 


29 


Kaz 102 


18 03 28.8 


67 38 10 


SI 


0.136 


4.62 


30 


KUV 18217+6419 


18 21 57.3 


64 20 36 


SI 


0.297 


4.05 


31 


MKN 885 


16 29 48.3 


67 22 42 


81.5 


0.026 


3.85 


32 


MKN 876 


16 13 57.2 


65 43 10 


SI 


0.129 


2.87 


33 


3C 390.3 


18 42 09.0 


79 46 17 


SI 


0.057 


4.24 


34 


MKN 290 


15 35 52.3 


57 54 09 


SI 


0.030 


1.72 


35 


MKN 279 


13 53 03.4 


69 18 30 


SI. 5 


0.0297 


1.78 


36 


H 1318+692 


13 20 24.6 


69 00 12 


SI 


0.068 


1.75 


37 


H 1419+480 


14 21 29.6 


47 47 27 


SI. 5 


0.072 


1.65 


38 


H 1320+551 


13 22 49.2 


54 55 29 


SI 


0.064 


1.36 


39 


PG 0804+761 


08 10 58.5 


76 02 43 


SI 


0.100 


2.97 


40 


MKN 506 


17 22 39.9 


30 52 53 


S1.5 


0.043 


3.26 


41 


MRK 6 


06 52 12.3 


74 25 37 


SI. 5 


0.0197 


6.39 


42 


H 1537+339'= 


15 39 52.2 


33 49 31 


SI 


0.329 


2.06 


43 


MKN 478 


14 42 07.5 


35 26 23 


SI 


0.077 


1.03 


44 


MKN 464 


13 55 53.5 


38 34 29 


SI. 5 


0.051 


0.96 


45 


NGC 5033 


13 13 27.0 


36 35 39 


LINER? 


0.0047 


1.03 


46 


PG 1425+267 


14 27 35.7 


26 32 14 


QSO 


0.366 


1.73 


47 


NGC 4151 


12 10 32.6 


39 24 21 


SI. 5 


0.0051 


1.98 


48 


AKN 564 


22 42 39.3 


29 43 31 


S1.9 


0.0247 


6.40 


49 


NGC 5548 


14 17 59.5 


25 08 12 


SI. 5 


0.0177 


1.69 



''SI : Scyfcrt 1, S2 : Seyfert 2 , SI. 5 : Seyfert 1.5, S1.9 : Seyfert 1.9 
"No data with ASCA or XMM-Newton observation 
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Table 3. Results - Sample 1 







Flux^ 




p 




Flux'^ 


Ohv ('rpfl'l 








HEAO-1 






I J-U IjLll J 


A<in A-XMM 






1 


IC 4329A 


8.2 


5.6 


-I 7. 1 +0.011 
— 012 


35.8+°-^ 


15.7 


XMM 


12.7/9.90 


2 


3C 273 


7.5 


540. 


1 fi^S+U-"29 

l.UOO_g 


0+9.7 


9.4 


XMM 


8.14/5.96 


3 


NGC 2992 


7.4 


1.0 


1.7 


90 


7.4 


GiOO 


59.2/27.1 


4 


NGC 5506 


5.5 


0.80 


1 721 

— 049 


323tr9 


12.0 


XMM 


13.1/9.92 


5 


NGC 526A 


5.2 


5.3 


1.61±0.02 


115±3 


3.5 


T97 


43.4/51.0 


6 


NGC 7582(*) 


3.9 


0.86 


1.95 


17440 


1.97 


TuOO 


56.4/52.2 


7 


ESQ 198-G024 


4.8 


23. 


1.77±0.03 


<3.2 


1.1 


Gu03 


0/6.8 


8 


MCG -6-30-15 


4.9 


0.73 


, Q2+0.02 


1 7+0.4 

^■'-0-3 


4.6 


R97 


147/ 


9 


MKN 509 


4.3 


13. 


1 404+0-029 

i.4y4_(j 028 


0+4-8 


3.4 


XMM 


24.7/0 


10 


3C 120 


4.9 


12. 


^•UU_o 02 


4.6t°o;^ 


4.4 


G98 


47.4/ 


11 


NGC 7172 


3.7 


0.97 




819+1^-2 


3.7 


T97 


14.9/15.6 


12 


NGC 3783 


4.2 


0.87 


1-601^2 




8.5 


BL02 


37.7/37.3 


13 


MKN 926 


4.0 


21.0 


l.Di^_Q 071 


0+26.3 


3.1 


XMM 


10.3/ 


14 


NGC 7469 


4.0 


2.3 


1 770+0-040 
^- ' "^-0-042 


0+7.3 


2.3 


XMM 


17.8/12.3 


15 


NGC 4593 


3.9 


0.86 


1 RQO+0-049 
l.D9^_Q 06 


0+8.8 


3.7 


XMM 


13.9/9.5 


16 


NGC 3227 


3.0 


0.25 


1.52 f 


6608 


0.8 


GO03(l) 


37.4/35.3 


17 


ESC 141-G55 


3.8 


12. 


1.72±0.06 


< 1 


2.7 


GO03(2) 


57.0/57.6 


18 


IIIZw2 


3.7 


76. 


1 7^+0.018 


0+0-21 


0.7 


XMM 


7.32/10.1 


19 


Fairall 49 


3.6 


3.5 


1 Q8q+0-044 
i.»oo_Q 044 


92.2+5-^ 


1.3 


ASCA 


20.2/34.5 


20 


3C 445(*) 


2.3 


28. 


1 ins+0-i57 

^■lU8_o.ll8 


1321111^ 


0.7 


XMM 


0/15.3 


21 


MKN 1152 


3.2 


21. 


1.672t»;°f, 


i3-8t^:^ 


0.6 


ASCA 


33.8/37.4 


22 


NGC 7314 


3.2 


0.16 


1.848lHli 


72tl 


4.0 


XMM 


42.3/30.4 


23 


MCG -1-24-12 


2.5 


3.1 


1 CQ+O.OT 
-i.oa_Q 03 




1.0 


Ma02 


76/36 


24 


ESQ 103-G35 


2.4 


2.5 


1 qr;c:+0.203 
i-yoo_0.193 


2028l:ig 


2.3 


XMM 


9.7/0 


25 


3A 0557-383(*) 


2.9 


10. 


1.8* 


37 


2.0 


T96 


40/40 


26 


MRK 590 


2.8 


4.5 


1.526«:?r 


0+12-3 


0.5 


XMM 


10.4/7.0 


27 


H 1846-786 


3.1 


44. 


1 Q2Q+0-025 
±.a^a_Q 024 


1 2+1-^ 


0.8 


ASCA 


59.9/76.9 


28 


Fairall 9 


2.9 


15. 


1 ye, +0.051 
l.(01_(j 044 


< 1 


1.2 


XMM 


28.7/27.7 



''Observed After absorption, (lO^n erg cm^^ s^l) 
''Intrinsic X-ray luminosity in 2 — 10 keV band (erg s~l) 
After subtracted with Galactic Njj 
''Reference — 

R97 : (Reynolds 1997), GiOO : (Gilli ct al. 2000), T97 : (Turner ct al. 1997), 

TuOO : (Turner et al. 2000), Gu03 : (Guainazzi 2003), G98 : (George et al. 1998), 

BL02 : (Blustin et al. 2002), Go03(l) : (Gondoin et al. 2003a), Go03(2) : (Gondoin et al. 2003b), 

Ma02 ; (Malizia et al. 2002), T96 ; (Turner et al. 1996), 

XMM : Our own data analysis with XMM-Newton observation, ASCA : Our own data analysis with ASCA observation 
"'MOS/PN of XMM-Newton, SIS/GIS of ASCA or MECS/PDS of Beppo-SAX observation 
*The source was excluded from the result of mean photon index. 

^Neutral dual absorber with covering fraction ~90%. please see Table. 2 in (Gondoin et al. 2003a). 

(*'For these sources, we have used the spectral fit results involving models with complicated absorption features. See Sect. 3.2 for 
details. 
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Table 4. Results - Sample 2 



Source ID 


Name 


Flux^ 


T h 


Photon Index 


Intrinsic Nu*^ 


Flux^ 


Obs.(ref)° 


Exp.time^ 






HEAO-1 


(lO'*'^ erg/s) 




(10^0 cm~2) 


ASCA-XMM 


(ks) 




29 


Kaz 102 (*) 


3.3 


232. 


1.9(fixed)f 


< 1 


.43 


ASCA 


16.8/19.7 


30 


KUV 18217+6419 


2.0 


594. 


1 906+°-°" 


Q+0.93 


1.6 


ASCA 


40.0/43.9 


31 


MKN 885 


2.3 


3.6 


1 nKo+0.069 


16.8lli 


0.6 


ASCA 


18.5/21.5 


32 


MKN 876 


0.97 


43. 


1 071 +0-035 
-^-0.034 


0+0.27 


0.49 


XMM 


7.2/2.6 


33 


3C 390.3 


1.1 


8.5 


1 fiqc;+().()22 
l.DOO_(, „21 


2 9+1-" 


1.8 


ASCA 


39.1/46.7 


34 


MKN 290 


1.5 


3.2 




4 4+2-0 


1.0 


ASCA 


40.1/45.5 


35 


MKN 279 


2.5 


5.0 


1. (D'l_g ggg 


0+4.4 


1.7 


XMM 


29.2/26.3 


36 


H 1318+692 


1.2 


13. 


1 000 + 0.208 


2.31^^3^ 


.15 


ASCA 


19.7/22.0 


37 


H 1419+480 


2.0 


26. 


1 QQ7+0.033 
■I-OJ' _0.031 




0.7 


XMM 


9.9/7.2 


38 


H 1320+551 


2.1 


22. 


1 679+0-254 
l.O(»_„ 203 




.10 


ASCA 


10.1/11.6 


39 


PG 0804+761 


1.1 


28. 


9 -187+0. 027 
-0.026 




1.0 


ASCA 


41.8/47.6 


40 


MKN 506 


0.91 


4.0 


-1 QQ9+0.072 




.63 


ASCA 


16.5/18.2 


41 


MRK 6(*) 


3.2 


3.9 


1.81 (fixed)* 


1043 


1.3 


XMM 


18.4/26.0 


42 


H 1537+3398 


3.6 


1350. 


1.8f 


< 1 


** 


Missing 




43 


MKN 478 


0.78 


11. 


2.213i0-f2 


1.6ll2« 


.27 


XMM 


26.2/24.8 


44 


MKN 464 


3.9 


26. 


1 ^Qn+0-215 
-i.utiu_g 208 


66.11^0 


.27 


XMM 


7.5/5.5 


45 


NGC 5033 


3.0 


0.15 


1.668lO;Oo^^ 




.26 


ASCA 


35.7/38.8 


46 


PG 1425+267 


1.1 


595. 






.28 


ASCA 


38.1/41.7 


47 


NGC 4151 


6.9 


0.53 


1.93 


340 


4.4 


Sc03 


110/91 


48 


AKN 564 


0.71 


0.98 


9 -1Q9+0.051 

Z.l»Z_g g>,j 


0+7.9 


1.5 


XMM 


11.8/0 


49 


NGC 5548 


3.0 


2.1 


1 ^65+0-016 


0+3.2 


5.3 


XMM 


27.3/19.8 



^Observed After absorption, (10^^^ erg cm^^ s^^) 

''Intrinsic X-ray luminosity in 2 — 10 IceV band (erg s"') 

After subtracted with Galactic Nh 

''Reference - 
Sc03 : (Schurch ct al. 2003), 

XMM : Our own data analysis with XMM-Newton observation, ASCA : Our own data analysis with ASCA observation 
"'MOS/PN of XMM-Newton, SIS/GIS of ASCA or MECS/PDS of Beppo-SAX observation 
*The source was excluded from the result of mean photon index. 
SNo data with ASCA or XMM-Ne-wton observation 

(*^For these sources, we have used the spectral fit results involving models with complicated absorption features. See Sect. 3.2 for details. 
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Table 5. Spectal Results for Soft excess 



Source ID 


Name 


ivioaei 


r or Temp. 


Dual absorber Njj 
and Covering fraction^ 
(10^0 cm-2) 


Intrinsic Nfj"^ 
(1020 cm-2) 


\ jV matrix 


1 


IC 4329A 


blackbody + 4 egde 


188+^ eV 






1704/1387 


2 


3C 273 


2 power-law 


3.768+_%T 






881/817 


4 


NGC 5506 


mekal + dual absorber 


86+1° 


220l (95%) 


94111 


780/747 


9 


MKN 509 


2 power-law 


o 412 + 0.121 






465/427 


13 


MKN 926 


2 power-law 








251/221 


14 


NGC 7469 


2 power-law 








655/607 


15 


NGC 4593 


bremss -|- edge 


358t^g eV 






376/303 


19 


Fairall 49 


broken power-law 


4 646+°-*^° 

4-040_Q 






1575/1415 


20 


3C 445(*) 


dual absorber 




132lt\lt (90%) 




202/159 


22 


NGC 7314 


blackbody 


55t_l cV 






1404/1151 


24 


ESO 103-G35 


blackbody -|- 2 dual absorber 


771^2 eV 


1145lg^ (94%) 


168li« 


185/157 










10321^83 (95%) 






26 


MRK 590 


2 power-law 






403/365 


28 


Fairall 9 


2 power-law 


Q QQ9+0.276 
o.oa^_Q 216 






881/789 


32 


MKN 876 


2 power-law 


^ 1 C9+0.353 
^■-'■"^-0.564 






202/156^ 


35 


MKN 279 


bremss + edge 








663/563 


37 


H 1419+480 


blackbody -1- edge 


180t?| eV 






700/587 


43 


MKN 478 


bremss 








321/275 


44 


MKN 464 


2 power-law -|- edge 


c f.r,7+0.486 
o.uzi _Q g26 






263/225 


48 


AKN 564 


2 power-law 


4Qjt:+0.211 
4:.yiO_Q 244 






286/234 


49 


NGC 5548 


blackbody -|- 3 edge 


187t J eV 






842/679 



*1. 5-3.0 keV band was extracted for spectral fitting. 
''Reflection fraction 
After subtracted with Galactic Nh 

(*'For these sources, we have used the spectral fit results involving models with complicated absorption features. See Sect. 3.2 for 
details. 
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Table 6. Results of Iron Line 



oource ID 


Name 


line energy 


normalization 


Equivalent width 






(keV) 


(MOS or CIS) 


(MOS or GIS) 


1 


IC 4329A 


6 382+?,-03i 
6.435+;!-"5?, 

— O.OoO 


1 46+"1?c-4 


76 eV 


4 


NGC 5506 


7.29+?-3«c-5 

— 5.60 


52 cV 


9 


MKN 509 


6.302+°-202 

0.206 


4.96l5g|e-5 


125 eV 


14 


NGC 7469 


6.458+°-2§! 

" 0.286 


2.96+^536-5 


103 eV 


15 


NGC 4593 


6 296+"'°'''^ 


5.96I2 87*^-5 


138 eV 


19 


Fairall 49 


6.497_„ 


2.69^o;?3c-5 


132 eV 


20 


3C 445{*) 


D.O(0_Q 240 


2.56+||;|c-5 


178 eV 


21 


MKN 1152 




l-37lg;gge-5 


202 eV 


22 


NGC 7314 


fi OQ1 +0.074 

D.091_Q g4g 


4.49t?;i^e-5 


98 eV 


26 


MRK 590 


6 ocn+0.060 


1.55j:f|40-5 


253 eV 


28 


Fairall 9 


6 434+"'i44 

"•^•J^-1.983 


9.69lfo^c-5 


605 eV 


31 


MKN 885 


6.056_,j 248 




418 cV 


33 


3C 390.3 




4.60li||e-5 


215 eV 


34 


MKN 290 


6.484l°:li 


3.13 2 o7^~^ 


273 eV 


35 


MKN 279 


b.480_o QgQ 


1.63li;26g_5 


86 eV 


37 


H 1419+480 


f. +0.029 


8.65l^^g\°e-6 


109 eV 


43 


MKN 478 


c CC1 +0.f>28 
fi-661-().4(i2 


1.32 Q "^^6-5 


506 cV 


48 


AKN 564 


D.4U0_Q 405 


1.7lt«5ie-4 


1400 eV 


49 


NGC 5548 


6 407+0-050 
D.4Ur_Q 04g 


2.82tJ;«2e-5 


26 eV 



(*)For these sources, we have used the spectral fit results involving models with 
complicated absorption features. See Sect. 3.2 for details. 



Table 7. The Best Fit Parameters of the HXLF Models 



Sample 




log L,^ 


71 


72 


All= 

Unabsorbed 
Absorbed 


3.68±0.55 
1.98±0.36 
1.60±0.41 


44 02+"-^** 
44.00{fix) 
44.00(fix) 


n q'3+0.20 

U.yj_o 29 
'J- '^-0.19 
^•^^-0.19 


9 CI +0.53 

^•■^^-0.22 
o 04+0.90 
•J-J'*_o gg 


U03 


2.64±0.18 


44.11±0.23 


0.93±0.13 


2.23±0.15 



Note. — Errors are 1 a for a single parameter. 

''In units of IQ-^h^Q Mpc'^. 

''In units of h^Q erg s~^. 

''only 2 < 0.2 and log Wh < 24 sources. 
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